He II Ly/3 Gunn-Peterson Absorption: New HST Observations, 

and Theoretical Expectations 

David Syphers 1 , Scott F. Anderson 2 , Wei Zheng 3 , Britton Smith 4 , Matthew Pieri 1 , Gerard 
A. Kriss 3 ' 5 , Avery Meiksin 6 , Donald P. Schneider 7 , J. Michael Shull 1 , Donald G. York 8,9 

ABSTRACT 

Observations of He II Lya Gunn-Peterson troughs have proved to be a valu- 
able probe of the epoch of helium reionization at z ~ 3. Since this optical depth 
can become unmeasurably large even for modest He II fractions, various alter- 
nate techniques have been proposed to push to higher redshift, and among the 
more promising is looking at higher order Lyman-series troughs. We here report 
four new observations of the He II Ly/3 trough, including new data on the only 
sightline with a prior Ly/3 observation. However, the effective optical depth ratio 
T e ff,/3 / V e ff iQ is not simply predicted by fp\pf f a X a = 0.16, and we analyze cosmo- 
logical simulations to find that the correct ratio for helium at z ~ 3 is ~0.35. In 
one case we infer r e ff )a > 8.8, strong evidence that helium was not fully reion- 
ized at z = 3.2-3.5, in agreement with previous measurements suggesting a later 
completion of reionization. 
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Introduction 



The full reionization of helium at z ~ 3 was a major step in the evolution of the 
intergalactic medium (IGM). He II requires four-Rydberg photons for ionization, and such 
hard photons are believed to be produced primarily in active galactic nuclei (AGN), rather 
than stars. Thus the helium reionizati on epoch was delayed versus that of hydrogen (and 
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20091 ; ISyphers et al.ll201ll ). Helium reionization was not only a major change for the IGM 



state, but it also injected substantial heat into the IGM, and thus affected hydrogen as well. 

Various indirect methods of constraining helium reionization via the H I Lya forest 
have been used, either by examining thermal broadening of Lyq forest featu r es associated 
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with redshift is subj ect to considerable controversy (IMcDonald et al.ll200ll ; lKim et al.ll2002al ; 
Meiksin et al.l 120101 ) . and while t here appears to be a real opacity dip at z ~ 3.2, its inter- 
pretation is not straightforward (IBolton et al.l 12009c iMcQuinn et al.ll2009l ). Measurements 
of IGM metal line widths could break the thermal/non-thermal broadening degeneracy and 
yield firmer indirect co nstraints, but these are as yet observationally impossible in the nec- 
essary density regime (IMeiksin et al.ll2010l ). Optical depth ratios comparing metal species 
whose ionization potentials straddle the He II Lyman limit (54.4 eV, 228 A) may show a 
change in the UV background in z ~ 3 -4 that is possibly associated with helium r eionization 



(|Songailalll998l: lAgafonova et al.l 120071 ) . but this observational result is disputed (IKim et al. 
2002bl ; lAguirre et al.l 12004 ). and such a change is not predicted in some models due to inho - 
mogeneities in the UV background or in metallicity ( iFurlanettol 120091 ; IBolton fc Vielll201ll ). 



The best constraints on helium reionization hav e come from direct H e II optical depth 
measurements, bot h at higher redshift (3.2 < z "> 3.9;ISvphers et al.ll2011l) and lower redshift 



(2.4 < z < 2.9; e.g- IShull et al. 



2010; 



Worseck et al.l201ll ). However, there are regions of the 



Lya Gunn-Peterson trough ( IGunn fc Petersonlll965l ) with no detectable flux transmission in 
some individual spectra and multispectra composites. Because the Gunn-Peterson optical 
depth of a transition depends on oscillator strength and wavelength, n oc /jAj, higher- 
order transitions have smaller opacities. As a result, He II fractions that are black in Lya 
(giving only lower limits on the optical depth) can have measurable non-zero flux in higher- 
order trou^hs : _rhis_method of looking at highe r-order troughs has been used with H I 
(e.g., iBecker et al.l 120011 ) and proposed for He II (jMcQuinnl 120091 ) . To date there has been 
only one He II Ly/3 optical depth measurement, for HE2347-4342 (henceforth HE2347; z = 
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2.9), a challenging e xtraction of noisy data on the short-wavelength SiC detectors of FUSE 
flzheng et al"]l2004bh . 

In this paper, we compare He II Ly/3 absorption, A^ = 256.317 A, to He II Lya, 
X a = 303.782 A, in thr ee new sightlines, with data from the Cosmic Origins Spectrograph 
(COS; lOsterman et al.ll201ll ). We also analyze the COS data on Ly/3 in the HE2347 sightline. 
In Section |2] we describe the new observations of He II quasars that allow Ly/3 trough analysis, 
and in Section [3] we discuss our optical depth measurements and the methodology used. 

It is important to note that the simple relationship we have given above for optical 
depths of various transitions is valid for effective optical dep ths only under the assu mption 
of a homogenous IGM, and does not hold in a clumpy IGM (jOh fc Fur lanettol 120051 ). While 
this inhomogeneity effect has been taken into account in recent hydrogen Gunn-Peterson 
studies, it has been overlooked thus far in the sparse helium reionization literature dealing 
with higher order transitions. In Section H] we consider the relationship between effective r a 
and ra in the context of analytical models and cosmological simulations. In Section [5] we 
discuss our new measurements in light of the opacity relationship established, and consider 
the implications for helium reionization. We conclude in Section El 



2. Observations 

In order to probe He II Ly/3, quasars must have sufficient FUV flux, and also lie in 
a specific redshift range. We need to observe the Ly/3 Gunn-Peterson trough at a redshift 
where it is not completely obscured by a black Lya Gunn-Peterson trough; in practice this 
means z a < 2.8, where the mean Lya optical depth has dropped to ~l-2. This constraint 
requires that the quasars are at zqso ^ 3.5. At the other end, the effective area of COS drops 
sharply for A < 1150 A, requiring the flux at the He II break f\ : break > 10~ 16 erg s _1 cm -2 A -1 
for quasars at z > 3.3, and substantially larger at lower redshifts. This is observationally 
challenging, as very few z ~ 3 quasars have any detectable FUV fl ux at all, although t he 



lists of such He II quasars have recently been rapidly expanding (jSyphers et al.l l2009al lb: 



Worseck et al.ll2011t Syphers et al., in prep). 



Our sample in this paper consists of four He II quasars observed with COS that fulfill 
these flux and redshift requirements. The first was observed in an ongoing campaign to verify 
more He II quasars, i/5Tprogram GO 12178 (Syphers et al., in prep.), with a relatively short 
reconnaissance spectrum: SDSSJ0915+4756 (henceforth SDSS0915, z = 3.34). Two more 
were observed in a followup program, HST GO 12249 (Zheng et al., in prep.), with long expo- 
sures: SDSSJ1253+6817 (henceforth SDSS1253, z = 3.47), and SDSSJ2346-0016 (henceforth 
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SDSS2346 , z = 3.51). SDSS09 15 is newly verifi ed as a He II quasa r, while SDSS1253 was 
verified in lSvphers et al1l2009al and SDSS2346 in lzheng et aDl2004al Below 1100 A the COS 
effective area is very small, so for quasars at lower redshifts (down to z ~ 2.7), the required 
flux to m easure Ly(3 Gunn - Peterson is /a, break ^ few xl0~ 15 erg s _1 cm -2 A -1 . Currently, 
HE2347 (IShull et al.l 120101 ) uniquely satisfies this requirement, and completes our present 
sample. 

We observed our new targets with COS/G140L (R ~ 2000-3000), and adopt 7 pixels 
(0.56 A) resolution element. This is approximately the FWHM of the line-spread 

function (LSF), although the L SF is non- G aussian with substantial wings, due primarily to 
mirror polishing errors on HST (lKrissll201ll ). The data were reduced using CALCOS 2.13.6, 
and we coadded individual exposures with custom software that flatfields the data and takes 
into account various detector defects. For further details on the GO 12178 observations 
and the removal of geocoronal emission lines, see Syphers et al. (in prep.), and on the 
12249 observations, see Zheng et al. (in prep.). We are interested in broad effective optical 
depth measurements with flux averaged over all available trough regions, but in any case, 
averaging over substantial optical depth fluctuations is unavoidable for G140L data, because 
single resolution elements have widths of ~100-160 km s" 1 . HE2347 was observed with 
COS/G130M (R ~ 16,000-20,000) and COS/G140L; for details on that observation, see 
Shull et al.l ( 2010 ). Since G130M currently covers only A > 1134 A, we use this data only 
for the Lya opacities, and use G140L for Ly/3. The observations are summarized in Tabled] 



HE2347 was previously observed with FUSE for 619 ks (IKriss et al.ll2001l ; IZheng et al. 

2004bh. COS/G140L h as ~10 cm 2 effective area at very short wavelengths (X < 1000 A; 



McCandliss et al.l 120101 ). which is similar to or slightly higher than FUSE (jSahnow et al. 
2000h. For HE23 47. these G140L data constitute only ~12 ks of the COS observations 



( IShull et al.ll2010l ). and there would therefore appear to be little hope of improving upon the 
FUSE measurement. However, the FUSE detector consisted of four different channels of two 
segments each, with two different coatings, and it often proved challenging to match flux 
calibration across segments, particularly for the lower-wavelength SiC channels. In addition, 
FUSE had a much larger stray and scattered-light background than COS, and although the 
dark current was about the same as COS per 2-D pixel, it was on average nearly 50 times 
larger than COS/G140L per Angstrom in the 1-D spectrum. We thus deem it worthwhile to 
analyze the COS HE2347 data to see what limits they set on t@. 

For the three new G140L observations, we fit the quasar continua with power laws over 
the entire available wavelength range, from the He II Lya break to a S/N cutoff in the 
red (typically ne ar 1900 A). We fi rst deredden the spectr a with a Fitzpatrick fc Massa UV 
extinction curve ( jFitzpatricklll999l ). using E(B — V) from ISchlegel et al.l (119981 ). For the fit, 
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we specify regions that appear to be largely free of H I Lya absorption to initialize, and 
then iteratively use all consistent continuum points. HE2347 has much higher-resolution 
G130M data, as well as much better S/N, which allows the fit regions to be defined quite 
precisely. The G130M data also allow a careful fit of Galactic H I Lya absorption, and thus 
an E (B — V). W e there fore use the power law index and normalization, and the extinction, 
from IShull et al.l (120101 ) for this object. The continuum fits are shown in Figure [TJ along 



with the quasar proximity zones, which we exclude from our optical depth calculations. 



3. Optical Depth Measurements 

Our goal of estimating He II Lya effective optical depths is reached in three steps: (1) we 
measure a raw He II Ly/3 r e fj directly from each spectrum, (2) we correct this measurement 
for contaminating lower-redshift He II Lyct r c g by one of three different methods, and (3) 
we predict r eff Q , using our corrected r e g^ and a conversion factor found from cosmological 
simulations. We describe the first two steps in this section, and the last step in section HJ 



3.1. Optical Depth Methodology 

Calculating optical depths in very low-count regimes, such as black troughs or low- 
sensitivity portions of the detector, requires care. We here set forth a methodology appro- 
priate for such calculations, which not only gives good results for individual spectra, but 
enables comparison between different spectra. 

The COS detector has a very low background that is dominated by dark current 
(~10~ 4 counts s~ l pix -1 , for a pixel in the 1-D spectrum). Nonetheless, in Gunn- Peterson 
troughs with near- zero flux, this background can be substantially larger than the signal. We 
are thus measuring signal counts, n s , and background counts, n&, in a regime where we can 
have both n s <C rib and n s ~ 0. The CALCOS pipeline takes n s = n Q bs — (nb) in each pixel, 
where the average background is found from regions offset in the cross-dispersion direction 
to unexposed portions of the detector, and smoothed over 100 pixels. Because the low back- 
ground means that often n obs < (n b ) (or indeed n obs = 0) in any given pixel, we often have 
the unphysical result that n s < 0. Thu s the nonparametric bootstrap median method we 



developed for use with ACS prism data ( jSyphers et al.l 120 111 ) is not useful here. 



Even assuming a simple parametrization (e.g., Poissonian) does not lead t o one clearly 



corre ct way of determining confidence intervals when near a physical boundary (IMandelkern 



20021 ). However, the most important thing is to choose a consistent method for analyzing all 
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relevant observations and simulated spectr a , and for this we adopt the "unified" frequentist 
confidence intervals of iFeldman fc Cousins! (119981 ). which use likelihood ratio ordering. This 
is a standard in the particle physics community when working near physical boundaries, 
and in our simulations of COS spectra, it generally has acceptable coverage even for very 
high optical depths. Feldman- Cousins confidence intervals also have the desirable property 
that they are "unified", i.e., as the signal is reduced compared to the background and zero 
signal becomes included in the lower limit, the intervals naturally transition from two-sided 
to one-sided intervals with the same coverage. For classical Neyman confidence intervals, if 
the decision whether or not to quote one-sided or two-sid ed intervals is based on the data 
(as it typically is in astronomy; IFeldman fc Cousins! Il998l call this "flip-flopping"), then the 
coverage from repeated experiments is not the claimed coverage. 

One important aspect to note is that because Feldman- Cousins confidence intervals are 
unified, they are not central (that is, repeated measurements yield confidence intervals which 
do not necessarily have equal probability of lying above or below the true parameter value). 
They do transition into central intervals for large signal counts, but we caution that when 
we deal with high optical depths, we are of course in the low-count regime. 

For troughs where the total signal count is formally negative, we recommend quoting 
two values in addit ion to the confidenc e inter val. First, the "sensitivity" of the experi- 
ment, as defined by IFeldman fc Cousins! (Il998l ). which is the average upper limit obtained 
for many experiments with the given Poissonian backgrou nd but zero signal. Se cond, simi- 
larly motivated, the "detection upper limit," as defined by lKashyap et al. As with 
the sensitivity, the detection upper limit depends purely on the detector, rather than the 
source, and is defined as the maximum intensity a source can have without having at least 
a probability (3 of being detected at a significance level a. Note that we use frequentist 
Feldman- Cousins c onfidence interv a ls for both calculations, rather than Bayesian credible 
intervals as used by iKashyap et al.l (12010! ) for the latter quantity. 



We estimate the expected number of counts, absent any absorption, by extrapolat- 
ing the quasar continuum fit into the Gunn-Peterson trough, and multiplying by the flat- 
field. The primary flat field features are shadows made by the grid wires above the detector 
(jOsterman et al.ll201ll ). which reduce the number of counts seen in those regions. Multi- 
plying by the flatfield puts these shadows into the extrapolated quasar continuum, similarly 
reducing the expected counts. We also take care to use the time-dependent sensitivity func- 



1 Note that Kashvap et al. ( 201dh use the terms "upper bound" and "upper limit" to distinguish between 
quantities that we prefer to call, more descriptively, "source upper limit" and "detection upper limit" . The 
latter is a property of the detector, and does not depend on the source examined. 
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tion to co nvert fluxes into count s for a given observation, since COS sensitivity does change 
over time (jOsterman et al.ll201ll ). We derive expected counts using the exact same detector 
regions as for the actual source counts and background for each individual exposure, which 
enables direct comparison between the two. Source and background counts from all relevant 
exposures are combined to derive a source counts confidence interval, which is converted to 
an optical depth confidence interval using the expected number of counts. 



3.2. Removing He II Ly« Contamination 



Our targets were chosen at redshifts high enough so that the Ly/3 trough is accessible 
to COS, but low enough so that it lies in the He II Lya forest, rather than the Lya Gunn- 
Peterson trough. Nonetheless, the opacity in the helium Lya forest is substantial at these 
redshifts (~l-2), and must be taken into account. We consider three different methods to 
do so, each with systematics that are large, but at least different from each other. 

Method 1 (Ml) is the technique commonly used when calculating H I Ly/3 Gunn- 
Peterson opacity, which is to remove the mean Lya opacity a t the given redshif t, using 



a fit calculated from spectra of other, lower-redshift quasars (e.g.. iBecker et al.ll200ll ). There 
have been such fits of mean He II opacity, but based only on a few sightlines (no more than 
two at any given redshift), and they a re thus subject to large systematic u ncertainty, as 
well as simply having very large scatter (jShull et al.ll2010l ; IWorseck et al.ll201ll ). In addition, 
the effective optical depths used in the fits have often been calculated for selected coherent 
regions of high or low r, rather than averaging over redshift bins without regard for chance 
structure, which distorts the results. Of course, the substantial variance between sightlines 
also implies that the average opacity, even perfectly determined, may not correspond well to 
the actual opa c ity in a specific sightline. We use the simple semianalytical r\ = 80 model of 
Worseck et al.l (120111 ) to estimate the average T e s,a(z), as it provides a good fit to the optical 
depths of three He II quasars in the Lya forest at z = 2.3-3.0. (This model is not a direct fit 
of the He II Lya data, relying as it does on a hydrogen distribution and r\ = Nu e u/ Nm, but 
it nonetheless fits the data quite well, and allows sensible extrapolation beyond where we 
have data — namely down to z = 2.1, where we need the opacity for the HE2347 sightline.) 

The accuracy of this model over small wavelength ranges in our data can be tested by 
comparing the M2 Lya measurements (free of Ly/3; see below) to model predictions over the 
same redshift range. The model underpredicts the opacity for the higher-redshift quasars, 
and overpredicts the opacity for HE2347, but the differences are relatively small compared to 
the opacity change seen in and out of the Ly/3 trough. For HE2347, SDSS0915, SDSS1253, 
and SDSS2346, respectively, the difference between the value measured just redward of the 
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Ly/3 trough and the model prediction for that region is At = —0.30, 0.30, 0.47, and 0.75. 
Such differences are of the magnitude expected given sightline variations in density and 
ionizing background (and thus 77), but demonstrate why this method is not perfect, even if 
the average opacity model were derived from more than a few quasars. 

Method 2 (M2) uses the opacity of the He II Lya forest just redward of the Ly/3 trough 
as an estimate for Lya contamination in the trough. This approach relies on the effective 
Lya optical depth not changing substantially over small redshift intervals (Az ~ 0.05-0.1), 
which is clearly not a perfect assumption. It has the advantage over the first method that we 
are comparing two regions of the same sightline, where we know reionization has completed 
(from an observational perspective, r < 5), as opposed to comparing different sightlines, 
which can complete reionization at somewhat different redshifts. On the other hand, it 
should tend to give a somewhat higher Lya opacity than the true value (because the opacity 
decreases with redshift on average), which would underestimate the Ly/3 opacity. In practice, 
we use the 20 A above the Ly/3 break to calculate the effective optical depth. 

Method 3 (M3) exploits a unique possibility when analyzing helium, not available to 
hydrogen Gunn-Peterson studies. Because we are working with FUV spectra in helium, we 
can use the optical spectrum to probe the same redshift in hydrogen, allowin g us to se e 



and a ccount for IGM structure. For HE2347, we use the VLT/UVES spectrum (IKim et al. 



200 7|), and for th e other quasars, we use spectra from the Sloan Digital Sky Survey (SDSS; 



York et al.l 120001 ) . This may appear to be the best method, since it is the only one to 
give information on the very same IGM region that is causing our He II Lya absorption. 
Nonetheless, it too has notable drawbacks in practice. It relies on converting H I opacity 
to He II, THeii/THi — Tj/4 (where the factor of 4 comes from the wavelength dependence). 
Typical values are r\ ~ 50-100, but the re is substantial var i ation in r± throughout the IGM, 



and this changes during reionization (jZheng et al.l l2004bl ; iFechner et al.l l2006t IShull et al. 



20101 ). As a result, this is our primary uncertainty in this method. This conversion becomes 
more complicated when the optical spectra do not resolve IGM structure, as is the case with 
our new targets and their SDSS spectra (R ~ 1800 in our region of interest), although not 
with the VLT spectrum of HE2347. 

We calibrate our effective r/ eS = 4 x (r eff H eii)/('7"eff, hi) by considering redshift regions of 
the He II Lya forest that are not Gunn-Peterson opaque, but are redward of the onset of Ly/3 
absorption (see Figure [2]). In these regions we can measure the effective Lya opacities with 
no Ly/3 contamination (in practice, we also avoid some regions of poor continuum fitting in 
the optical data near the O VI emission line). In all three sightlines for which we have SDSS 
data, this comparison gives similar r] e g, and we thus adopt the single value r) e s — 44 ± 5 for 
all our sightlines. This approach neglects the redshift evolution of rj, as well as continuum 
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uncertainties, which are hard to characterize for the SDSS data. There may be a weak trend 
of increasing r) e g with increasing redshift, but we caution that we are not measuring rj, which 
is defined point by point rather than averaged over a large region like 7/ e ff. So although our 
SDSS rj e g is consistent with , e.g., the mean (77) ~ 50 found at lower redshift in the HE2347 



sightline (IShull et al.ll2010l ). this probably has little significance. Indeed, for HE2347 itself 



we prefer 77^ = 26 ± 5, and use this value with the VLT data. 

An additional source of uncertainty in this method is continuum fitting the SDSS data, 
for which we use a well-established spline-fitting routine. We tested this routine on mock 
SDSS spectra of comparable S/N to our actual quasars, and found that it introduced errors 
on average ~10% (ranging from 1% to 20%) in the H I Lya r fj measurement, typically 
underestimating the opacity. This is a non-trivial source of error, but the systematic un- 
certainties in r) e ff dominate the error budget. The presence of H I Ly/3 and O VI quasar 
emission do complicate the optical spectrum, but the redshift region of the Ly/3 trough is 
largely unaffected, due to our exclusion of any proximity zone (see Figured]). 

Before converting r eff H i into r effj Heii, there is a small correction we take into account. 
Because this region of the H I Lya forest also contains higher-redshift H I Ly/3 lines as well, 
we remove that source of opacity. To do this, we measure the relevant H I Lya forest at 
higher redshift, and estimate the concomitant H I Ly/3 forest, using mean flux decrements 
and the effective H I opacity ratio r c s t p/r c s ta = 0.29 (as found in Section @J. The effect of 
this is shown in Figure [2J It might initially appear interesting to examine the He II Ly/3 
opacity without this last step, because using all H I opacity (including Ly/3) when converting 
to helium opacity should give a He II Ly/3 optical depth from only the diffuse IGM not 
detected in H I absorption. This also being opaque would be very strong evidence for Gunn- 
Peterson absorption, but in practice, the uncertainty in r/ eff makes such a test impractical. 



3.3. Optical Depth Results 

For each object, we exclude any possible proximity zone (as seen in He II Lya; see 
Figured]), where radiation from the quasar itself noticeably affects local optical depths. We 
then measure the opacity of the Ly/3 troughs (Figure [3]), using the methodology of Section 

The size of the Ly/3 trough we can measure is limited by two factors. For HE2347 and 
SDSS2346, we use the full trough from the onset of Ly/3 absorption down to the onset of 
Ly 7 absorption (243.027 A rest frame). For SDSS1253 and SDSS0915, our lower limits are 
determined by S/N, itself set by the rapidly declining response of the COS detector towards 
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shorter wavelengths. (We can test the difference of these two cutoffs in SDSS1253, where 
we have data of quality that is usable, albeit low, down to the Ljj break. Using this larger 
definition of the trough yields results very consistent with the S/N cut, although with larger 
error bars due to inclusion of very low-S/N data.) 

We remove geo coronal emission lines by using only data taken during orbital night in 
those regions affected (for details, see Syphers et al., in prep.) Despite this removal, a 
remaining concern might be that the very high flux of the geocoronal lines (particularly 
Lya) is large enough that scattered light from the lines will contaminate the spectrum far 
from the line wavelengths. To test this possibility, we calculate optical depths in troughs 
two ways, first using all data taken, and then using night-only data. (Although geocoronal 
Lya is present at all times, it is substantially weaker at night.) In every case the results are 
consistent, including for the highest optical depths, and thus scattered geocoronal emission 
is negligible. 

Results of our measurements are presented in Table [2j The columns of this table are as 
follows. Column (1) — target name. Column (2) — target redshift, taken from a combination 
of the onset of IGM He II Lya absorption and rest-frame UV emission lines seen in the optical, 
except for HE2347, for which the redshift is taken fro m [O III] A5007 (R. Simcoe, personal 



communication), and which agrees with O I A1302 (IReimers et al.l 119971 ). Column (3)— 
the redshift range used for the He II Ly/3 trough (Fig. [3]), avoiding any quasar proximity 
zone, the Ly7 trough, and extremely low-S/N areas. Column (4) — the measured He II Ly/3 
trough optical depth, with no correction for contaminating Lya opacity. Quoted errors 
are 68% Feldman- Cousins confidence intervals. Columns (5-7) — estimated optical depth for 
He II Lya in this redshift region, using methods Ml, M2, and M3, described in Section l3~2l 
Errors in column 7 are 68% random only. 



4. Predicted Optical Depth Ratios 



It is widely understood that (r) = (— ln(F/F n )) doe s not equal r e ff = — ln((F/F )) 
in an inhomogeneous IGM (e.g., iDixon fc Furlanettoll2009l ). but less widely appreciated is 
the implication that T a /rp ^ T c fj ia /r e fj i( g. As a result, when using T e g } p to estimate T c s jCe , we 
cannot simply use the ratio of oscillator strength and wavelength (fiXi), as has been done 
in the h eretofore sparse He II li terature on this topic, as well as early H I Gunn-Peterson 
studies. lOh fc Furlanettol (120051 ) pointed out this issue, and derived a ratio applicable for 
z ~ 6 hydrogen, but we need to reassess this for our z ~ 3 helium troughs. 



To characterize this effective optical depth ratio, we turn to the cosmological simulations 
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of ISmith et al.l ( 1201 ll ). These simulations use the Eulerian hydrodynamics + iV-body code 
Enzo, with adaptive mesh refinement turned off to allow uniformly good resolution in the 
low-density IGM. The simulation box was 50 h~ x comoving Mpc on a side, with 1024 3 
cells, a dark matter particle mas s of 7 x 10 6 h~ x M , and distributed feedback (simulation 
50_1024_2 of ISmith et al.l 1201 lh . Radiation post-processing for metal species was done 
with the UV background models of lHaardt fc Madaul (120011 ). assuming equilibrium; non- 
equilibrium ionization states of H and He were calculated as the simulati on ran. For details o n 
the feedback method, the ionization calculation, and other aspects, see ISmith et al.l ( 120111 ). 



We take 500 rays through this volume at each of several redshifts from z = 2.5-3.8, 
and create a high-resolution (R > 50,000) simulated spectrum for each ray. We verify that 
point-by-point, (r / 3(A)/r ct (A)) = 0.160 and (r 7 (A)/r a (A)) = 0.0557, showing that we have 
enough resolution to recover the fi\ ratio we expect in this case (Table [3]). 

To test how these spectra re produce the z ~ 3 forest , we co mpare the hydrogen r eff Q ,(z) 
to the observational values from IFaucher-Giguere et al.l (120081 ) . The simulations track the 
effective optical depth quite well, with the exception of the highest-redshift simulation point 
used (z ~ 3.7), which differs by 0.3 dex from the observed value, and a ~2cr discrepancy 
near the z ~ 3.2-3.3 opacity dip. This apparent observed dip (a devi ation from a power law 



fit) is not well understoo d, nor well reproduced by simulations (e.g., IFaucher-Giguere et al. 



2008 



Bolton et al.l 120091 ) . However, the opacity ratio of interest to us is fairly insensitive 
to minor disagreements such as this; over the redshift range z = 2.5-3.7, our hydrogen and 
helium optical depths change by a factor of ~4, but the ratio Tp/r a varies by only 7% for 
helium. We also compar ed the flux proba bility distribution function of the simulation with 
the observed values from lKim et al.l ( 120071 ). at z = 3.0. There is good agreement here except 
in the two highest flux bins (transmission T > 0.925), where some substantial discrepancy 
i s predicted based on proble ms finding the continuum observationally at higher redshift 



( IFaucher-Giguere et al.l 120081 ). 



The effective optical depth ratios predicted by the simulation, contrasted with the naive 
expectation, are shown in Tabled the colum ns of which are as follow s. Column 1 — transition. 
Column 2 — wavelength of the transi tion dRalchenko et al.l 120081 ). Column 3 — oscillator 



strength calculated from equation 8 of iMeiksinl ( 120091 ) (which agrees with lVerner et al.lll996 
to the precision that the latter reference quotes). Column 4 — the predicted optical depth 
ratio given by the ratio of /jAj. Column 5 — the ratio actually seen in simulations, for redshift 
z ~ 3.3. 



The optical depth ratios do evolve with redshift, albeit weakly, trending towards the 
point-by-point value as the mean opacity drops. We find for He II at redshifts z = 2.4-3.8, 



-12- 



to within 2% for Ly/3 and 3% for Ly7, 




0.190 x (1 + z) 



0.434 



(la) 



Q 



0.0663 X (1 + z) 



0.707 



(lb) 



a — 



Assuming a constant ratio over this redshift range is not a bad approximation for Ly/3, 
with deviations of <7% from the specific simulation values, but there is some systematic 
bias in doing this. We therefore use the redshift-dependent version, equation [Taj in our 
helium r e ff jQ estimates in Table HJ The Ly7 ratio is constant over this redshift range to 
within 12%. Large-scale structure variations cause a standard deviation over all sightlines 
of ~8-10% in the Ly/3 ratio, and ~12-19% in the Ly7 ratio, so the ratios could be ~10% 
wrong for a specific sightline. (Although Ly7 is not used in the current work, it is included 
for completeness, as it may be accessible to future observation. In particular, new modes of 
COS with higher blue throughput could enable this analysis.) 

These results used noise-free, normalized simulated spectra, so to determine if noise 
has an impact on these results, we create versions that are realistic mock COS spectra. 
We put in a power law continuum, redden the spectrum, convolve with the COS LSF, find 
the expected counts in each pixel, and generate a random realization of each spectrum 
assuming Poissonian statistics. We find that noise from random sampling has a negligible 
effect on our resultant ratio, as one would expect for unbiased noise and a sufficiently large 
sample. However, low fluxes and short observing times can affect the observed ratio; e.g., 
for /1200 A = 5 x 10~ 17 erg s _1 cm -2 A -1 and an exposure time of 3 ks, the mock observation 
ratio differs from the true simulation ratio by ~20%. However, for all the quasars in this 
paper, the exposure times and fluxes are sufficient to keep this bias <2%. 

While the primary purpose of the simulations is calculating the z ~ 3 He II T e ^^/T e Q^ a 
ratio, we also use them to find that for H I at z ~ 0, T e s^/T e s >a — 0.29. We need the latter 
quantity to remove low- redshift H I Ly/3 contamination from our H I Lya T e g measurement, 
as detailed in Section 13.21 



We apply the effective optical depth ratios to the Ly/3 measurements to predict he- 
lium T e s >a , correcting the Ly/3 optical depths using each of the three methods discussed in 
section 13.21 The results are presented in Table HI the columns of which are as follows. Col- 
umn (1) — target name. Columns (2,4,6) — He II Ly/3 optical depth measurements, taking 
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the raw optical depths from Table [2] and correcting by methods Ml, M2, and M3, respec- 
tively. Columns (3,5,7) — predicted He II Lya optical depths, using the Ly/3 measurements 
of Columns 2, 4, and 6, and the theoretical ratios derived in Section |H Errors quoted are 
68% random only Column 8 — measurement of the Lya optical depth in the same redshift 
range as the Ly/3 measurement. For the SDSS quasars, only 68% lower limits are quoted, as 
discussed below. 

In some cases, the predicted Lya optical depth is small enough that we should be able 
to measure it, allowing us to observationally find r cffi/3 /r cff a . We do this for HE2347, but 
unfortunately, the SDSS quasars are all contaminated by geocoronal O I A1302 in the Lya 
region corresponding to our Ly/3 measurements. Although O I is dramatically weaker at night 
(Figured] uses night-only data in that region), it still is strong enough to noticeably bias any 
measurements in a Gunn-Peterson trough. Table H] therefore quotes only 68% lower limits 
on the Lya optical depths for the SDSS quasars, because of this residual geocoronal flux. 
The HE2347 measurement and the lower limits on SDSS0915 and SDSS1253 are sufficient to 
show that M3 appears to be a very poor correction method, substantially underpredicting 
Lya opacity. This is likely due to the large uncertainty in translating r eff H i into Tef^Heii- 
Even if the simulations were wrong, and the effective optical depth ratio were as low as 
possible (0.16), the SDSS1253 data strongly disfavor M3. The measured Lya values are 
consistent with the predicted values for Ml and M2, using the simulation ratio. Given the 
systematic problems we have seen that Ml has with redshift, we regard M2 as perhaps the 
best estimate. This is supported by the good agreement between the M2 prediction and the 
actual measured Lya opacity in HE2347 (Table H]). In addition, the predictions of M2 are 
more conservative than those of Ml, and therefore, of our three methods, its estimates are 
the most conservative that fit the data. 

Using the M2 predictions, we set a lower limit on the Lya opacity of Tgff^eii > 8.8 in 
the SDSS2346 sightline at (z) = 3.38, and estimate r eff H eii = 4.361};^ in the SDSS1253 
sightline at (z) = 3.35 and r cff>H eii = 3.89±|^ in the SDSS0915 sightline at (z) = 3.29. For 
HE2347, at (z) = 2.78, we estimate r effiH eii = 2.51+°;49. 

Our results show that sightline variance, such as that between SDSS1253 and SDSS2346, 
is quite strong at high redshift, as it was already known to be at low redshift. It is diffi- 
cult to turn these mean optical depths into detailed information on the progress of helium 
reionization, however, in part because of the still-small number of sightlines, but also be 



cause theoretical m odels differ substantially at these redshifts (IDixon fc Furlanettol 12009 



Worseck et al.ll201ll ). Our own simulations were not intended for studying the evolution of 



helium reionization in detail. The large sightline variance itself m ay imply that helium reion- 



ization is incomplete at these redshifts (IFurlanetto fc Dixonll2010l ). but along some sightlines 



-14- 



it was already k nown to be incomplete at even lower redshifts (e.g., at z — 2.8 in HE2347; 



Shull et al.ll2010f ). Of course, the large majority of helium ionization and heat injection occur 



prior to the observational completion of reionization, which can be thought of as the end of 
black (r > 5) Gunn-Peterson troughs. 

The He II optical depth expectations for a uniform IGM (equation [2j below) have led 
to pessimism about being able to observe helium reionization by this method, but these 
expec ted values are m uch higher than they would be in a more realistic inhomogeneous 
IGM. iMcQuinnl ( 120091 ) discusses some aspects of this , and our simul ations show it well in 



our favore d metric, r e ff. Fo r the WMAP7 parameters (IKomatsu et al 



2011) and the helium 



fraction of ISteigmanl ( 120071 ). one obtains for a uniform IGM (jGunn &: Peterson! Il965f ) 



T"GP,HeII = 3.49 



uo-v V i 

1.0209 



Y P 



0.2486 / V 0.0456 



1 + 0.0209- (4/(1 + z)f 



h 



0.704 / V 0.272 



-1/2 



l + z 



3/2 



(2) 



where the last term is included to account for the cosmological constant (a small but non- 
negligible correction at z = 2-3), A b = Pb/(pb) is the baryon overdensity, and Y p is the 
helium mass fraction. 

For a uniform IGM, therefore, the helium Gunn-Peterson trough becomes saturated to 
the point that nonzero flux can no longer be measured at XHeii ~ 2 x 10~ 3 , even for a bright 



( Furlanetto 


2009; 


McQuinn 


2009) 



simulation, while not intended to match the details of helium reionization, provides a clear 
quantitative example of the effects of density variation. For example, at an ionization fraction 
^Heii — 2.5 x 10~ 3 (the mass- and volume-weighted ionization fractions x^ eU ~ ig eII in this 
case), we find helium r e s t0l = 1.65, rather than the 8.7 predicted at mean density by equation 
[2j The lower value arises because in an inhomogeneous IGM, transmission is dominated by 
under dense regions. 



6. Conclusion 

We present He II Ly/3 Gunn-Peterson trough measurements for three new quasars, 
and analyze new data on the one quasar for which Ly/3 had been measured before. The Ly/3 
trough has lower opacity than Lya, and thus allows measurements at otherwise unobservable 
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He II fractions. We predict large but variable Lya optical depths at z ~ 3.2-3.5, from 
r cS ,a ~ 4 to >9. 

However, the lower opacity of Ly/3 is not simply predicted by combining the Lya opac- 
ity with the fi\i ratio, when dealing with effective optical depths. We use cosmological 
simulations to derive r^p/r^a ~ 0.35 for helium at z ~ 3. Where both Ly/3 and Lya 
measurements are possible, they suggest good agreement between our predictions and the 
observations, although systematic errors make predicting Lya from Ly/3 somewhat uncertain. 
The He II Ly7 trough, although a factor of three stronger in r e g than in a point-by-point 
measurement, will be challenging to use. The short wavelength span of this trough and its 
lower wavelength position (at lower detector sensitivity) will make the opacity measurements 
quite difficult, and the corrections will need to include both Lya and Ly/3, and thus be that 
much more uncertain. 

We present a trough optical depth calculation method suitable for COS (with confi- 
dence intervals appropriate for any photon-counting device). Because of the unified nature 
of Feldman- Cousins confidence intervals, they allow comparison between all sets of data us- 
ing this method, regardless of whether the intervals are one- or two-sided. This aspect is 
important for any theoretical studies that wish to compare helium reionization simulations 
to data from different surveys. 

He II Ly/3 Gunn- Peterson troughs appear to be an effective way to measure high optical 
depth regions, and new short- wavelength modes of COS (HST/GO programs 12501 and 
12505) will be able to make better measurements of this along several sightlines. 

We thank Tae-Sun Kim for providing the continuum-normalized VLT spectrum of 
HE2347. 
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Fig. 1. — Regions near He II Lya for the four quasars studied in this paper. The vertical 
dashed line indicates the onset of He II Lya absorption according to the quasar redshift, and 
the shaded region is the ionized proximity zone, identified in Lya. The redshifts covered by 
the proximity zone are excluded from Ly/3 analysis. Continuum fits are overplotted (solid 
lines), and use data from the He II Lya break to a S/N cutoff in the red (typically near 
1900 A). Residual geocoronal O I emission is visible near 1302 A (marked). 
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Fig. 2. — Continuum-normalized spectrum of SDSS2346, as an example of how the M3 
estimator is derived. The upper curve is the optical data for H I, and the lower curve (shown 
only for z > 2.805) is the UV data for He II. Solid horizontal lines indicate fluxes of 
and 1, while short-dashed horizontal lines indicate the mean measured flux for the given 
spectrum in the relevant shaded region. n e s is found by comparing these mean fluxes in the 
calibration region (the right shaded region; green in the online version), where there is only 
Lya absorption. The mean flux in the region where M3 is measured (the left shaded region; 
blue in the online version) is contaminated by H I Ly/3 absorption from higher redshift; when 
we correct for this, the mean flux moves up to the long-dash line. 
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Fig. 3. — He II Ly/3 regions of the four quasars used in this paper. Extrapolated continuum 
fits are overplotted (solid lines). The vertical dashed lines indicate the onset of He II Ly/3 or 
Ly7 absorption, and the shaded region indicates the Ly/3 trough as extracted. This trough 
avoids the proximity region identified in Figure [TJ The blue edge of the trough is set by the 
onset of the Ly7 trough in HE2347 and SDSS2346, and by very low-S/N regions in SDSS0915 
and SDSS1253. Ly/3 optical depths for the latter two are robust under ~10 A changes in 
the blue edge, because the sensitivity is so low there. 
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Table 1. Observations 



Target 


Grating a 


Obs. Date 


Exp. Time 

(s) 


HE2347-4342 


G130M/G140L 


2009 


Nov 5 


28,458/11,558 


SDSSJ0915+4756 


G140L 


2010 


Oct 6 


5,521 


SDSSJ1253+6817 


G140L 


2011 


May 5 


14,096 


SDSSJ2346-0016 


G140L 


2010 


Nov 29, Dec 4 


20,737 



a Of HST/COS. G130M covers A = 1134-1474 A, while G140L covers (at 
greatly varying sensitivity) A ~ 920-2000 A, where the lower limit is from 
the Galactic Lyman limit. Both gratings have multiple central wavelength 
settings that shift coverage, and the values quoted are extremes. 



Table 2. Measured He II Ly/3 Optical Depths 



Target 


Quasar Redshift 


He II Ly/3 redshift range 


Raw r cff>/3 


Ml T eff , Q 


M2 r effiQ 


M3 r effiQ 


HE2347-4342 
SDSSJ0915+4756 
SDSSJ1253+6817 
SDSSJ2346-0016 


2.887 
3.343 
3.470 
3.510 


2.684-2.880 
3.253-3.329 
3.292-3.414 
3.276-3.493 


i 7^+0.14 
L - '°-0.12 

3 46+ 1 - 37 
o qq+o.66 

[6.02, oo) a 


1.07 
1.66 
1.74 
1.79 


92 +o n 
2 10+ - 12 

9 /IQQ+0.084 
z.<±oo_ _ 080 

9 81 A+0.070 
z.O±^±_ 068 


1.374 ±0.044 
2.54 ±0.40 
3.66 ±0.52 
3.24 ±0.34 



Note. - - Method 1 (Ml) uses an average semianalytical approximation to r cff)Q (z), method 2 uses r eff Q from just 
above the Ly(3 break, and method 3 uses the H I Lya optical depth to predict r e fj jQ . See Section I3~2l for further details. 

a No finite best estimate exists for this quantity. The sensitivity is r s = 5.59 (68%), and the detector lower limit is 
T d n = 5.15 (a = 0.32, /3 min = 0.9). (See Section [37T1 of the text for discussion of sensitivity and detector lower limit.) 
Using the night spectrum only, we obtain t c q^ = 6.22^0^. 
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Table 3. Simulation Predictions for the Optical Depth Ratio 



Transition 


Wavelength 


Oscillator Strength (/) 


n{\)/r a {X) 


Teff,i/VcfT,a a 


Lya 


303.782 


0.4162 


1 


1 


Ly/3 


256.317 


0.07910 


0.160 


0.35 


Ly7 


243.027 


0.02899 


0.0557 


0.18 



a As found in simulations for z = 3.3; see Section HI These do not vary 
strongly with redshift. 



Table 4. Predicted He II Ly« Optical Depths 



Target 


T e fi,p, corr. Ml 


T e S,a, pred. Ml 


T e ff,/3, corr. M2 


Toff, a , pred. M2 


r e ff,/3, corr. M3 


T e s,a, pred. M3 


T e ff,a, measured 51 


HE2347-4342 
SDSSJ0915+4756 
SDSSJ1253+6817 
SDSSJ2346-0016 


0.68^| 

1- 8 0-0.69 
9 9 r+0.66 
z - zo -0.40 

[4.23, oo) 


9 n fi +0.43 
z - UD -0.37 

5.i4i?; 9 97 

6-25±i;S 
[11.59,oo) 


o.83i8:S 

1 5fi +0 - 67 
i.oo_ 41 

[3.16, 00) 


9 ci +0.55 
Z - Oi -0.49 

3.89l|;g 
4-36±i:fI 
[8.78, 00) 


o-38iS:iS 

0.92t ;« 

0.33tg£ 
[2.54,00) 


2.63^ 
[7.06, 00) 


2.393 ±0.015 

> 3.84 

> 3.55 

> 4.49 



Note. - - Ly/3 optical depths are corrected using the estimates of Table [2] and Lya optical depths are predicted using the ratio of Section |4j 
Errors quoted are random only. 

a Only lower limits exist for the three SDSS quasars, because of contaminating geocoronal O I emission. 



